Order-Degree Curves for Hypergeometric 
Creative Telescoping 
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ABSTRACT 

Creative telescoping applied to a bivariate proper hypergeo- 
metric term produces linear recurrence operators with poly- 
nomial coefficients, called telescopers. We provide bounds 
for the degrees of the polynomials appearing in these opera- 
tors. Our bounds are expressed as curves in the (r, d)-plane 
which assign to every order r a bound on the degree d of 
the telescopers. These curves are hyperbolas, which reflect 
the phenomenon that higher order telescopers tend to have 
lower degree, and vice versa. 

Categories and Subject Descriptors 

1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms 

General Terms 

Algorithms 

Keywords 

Symbolic Summation, Creative Telescoping, Degree Bounds 

1. INTRODUCTION 

We consider the problem of finding linear recurrence equa- 
tions with polynomial coefficients satisfied by a given defi- 
nite single sum over a proper hypergeometric term in two 
variables. This is one of the classical problems in symbolic 
summation. Zeilberger [16] showed that such a recurrence 
always exists, and proposed the algorithm now named after 
him for computing one [15, 17]. Also explicit bounds are 
known for the order of the recurrence satisfied by a given 
sum [14, 10, 4]. Little is known however about the degrees 
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of the polynomials appearing in the recurrence. These are 
investigated in the present paper. 

Ideally, we would like to be able to determine for a given 
sum all the pairs (r, d) such that the sum satisfies a linear 
recurrence of order r with polynomial coefficients of degree 
at most d. This is a hard question which we do not expect to 
have a simple answer. The results given below can be viewed 
as answers to simplified variants of the problem. One sim- 
plification is that we restrict the attention to the recurrences 
found by creative telescoping [17], called "telescopers" in the 
symbolic summation community (see Section 2 below for a 
definition) . The second simplification is that instead of try- 
ing to characterize all the pairs (r, d) , we confine ourselves 
with sufficient conditions. 

Our main results are thus formulas which provide bounds 
on the degree d of the polynomial coefficients in a telescoper, 
depending on its order r. The formulas describe curves in 
the (r, d)-plane with the property that for every integer point 
(r, d) above the curve, there is a telescoper of order r with 
polynomial coefficients of degree at most d. As the curves 
are hyperbolas, they reflect the phenomenon that higher or- 
der recurrence equations may have lower degree coefficients. 
This feature can be used to derive a complexity estimate 
according to which, at least in theory, computing the mini- 
mum order recurrence is more expensive than computing a 
recurrence with slightly higher order (but drastically smaller 
polynomial coefficients). This phenomenon is analogous to 
the situation in the differential case, which was first analyzed 
by Bostan et al. [5] for algebraic functions, and recently for 
integrals of hyperexponential terms by the authors [6] . 

Our analysis for non-rational proper hypergeometric input 
(Section 3) follows closely our analysis for the differential 
case [6]. It turns out that the summation case considered 
here is slightly easier than the differential case in that it re- 
quires fewer cases to distinguish and in that the resulting 
degree estimation formula is much simpler than its differ- 
ential analogue. For rational input (Section 4), we derive a 
degree estimation formula following Le's algorithm for com- 
puting telescopers of rational functions [2, 9]. 
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2. PROPER HYPERGEOMETRIC TERMS 
AND CREATIVE TELESCOPING 

Let I be a field of characteristic zero and let K(n, k) 
be the field of rational functions in n and k. We will be 
considering extension fields E of K(n, k) on which two iso- 
morphisms S„ and Sk are defined which commute with each 



other, leave every element of K fixed, and act on n and 
k via S„(n) = n + 1, S k {n) = n, S„(k) = k, S k (k) = 
k + 1. A hypergeometric term is an element /i of such an 
extension field E with S„(h)/h G K(n, fc) and Sk(h)/h G 
K(n, fc). Proper hypergeometric terms are hypergeometric 
terms which can be written in the form 

7 n k TT F(& m 77. ~t~ &mk + <2 m )F(f) m ?2 6 m 

m / 

PX ^ ^ T(w m n + iCA: + u!^)T(v m n - v' m k + v m ) ' 

(1) 

where p G K[n, fc], a;, y G K, M G N is fixed, a m , a' m , b m , 
b' m , Um, u' m , v m , v' m are nonnegative integers, a' m , b' m , u'm, 
uJJ, G K and the expressions x n , y k , and F(. . .) refer to 
elements of E on which S n and Sk act as suggested by the 
notation, e.g., 

S n (T(2n -k + 1)) = (2n-&+l)(2n-fc+2)r(2n - k + 1), 

S t; (r(2n - k + 1)) = ^^: r (2™ -k + 1). 

Throughout the paper the symbols p, x, y, M, a m , a' m , 
a m , . . . will be used with the meaning they have in (1). For 
fitting long formulas into the narrow columns of this layout, 
we also use the abbreviations 

A m := a m n H- c^ink + G* m7 B m := b m n + b m k -}- b m , 

Um := iimn + u' m k + u'm, Vm := v m n + v' m k + v' m . 

We follow the paradigm of creative telescoping. For a 
given hypergeometric term h as above, we want to determine 
polynomials Iq, . . . ,l r G K[n] (free of k, not all zero), and a 
rational function C G K(n, k) (possibly involving k, possibly 
zero), such that 

£ h + tiS n {h) + ■■■+ l r S r n {h) = S k (Ch) - Ch. 

In this case, the operator L := £o + £iS n + • • • + £rS„ G 
K[n][S„] is called a telescoper for h, and the rational func- 
tion C G K(n, k) is called a certificate for L (and h). The 
number r is called the order of L, and <f := maxj_ deg n ii is 
called its degree. If ft represents an actual sequence f(n, k), 
then a recurrence for the definite sum X^=o /( n > ^) can ^ e 
obtained from such a pair (L, C) as explained in the litera- 
ture on symbolic summation [12]. We shall not embark on 
the technical subtleties of this correspondence here but re- 
strict ourselves to analyzing of the set of all pairs (r, d) for 
which there exists a telescoper of order r and degree d. 
The following notation will be used. 

• For p G K[n, k] and m G N, let 

p m := p(p + l)(p + 2) • • • + m - 1) 
with the conventions p° = 1 and p 1 = p. 

• For p G K[n, fc], deg n p and deg fe p denote the degree 
of p with respect to n or k, respectively, degp without 
any subscript denotes the total degree of p. 

• For z G R, let z + := max{0, z}. 
With this notation, we have 



S n {h) S n (p) 

— ; = X 



n si-m D m 

m — l 



S k (h) _ Sfc(p) ^ Afr(V m - v' m y 



3. THE NON-RATIONAL CASE 

We consider in this section the case where h cannot be 
split into h = qho for q G K(n, k) and another hypergeomet- 
ric term ho with Sk(ho)/ho = 1. Informally, this means that 
we exclude terms h where y — 1 and every F-term involving 
k can be cancelled against another one to some rational func- 
tion. Those terms are treated separately in Section 4 below. 
If h cannot be split as indicated, then also Ch cannot be 
split in this way, for any rational function C G K(n, k). In 
particular, we can then not have Ch G K and therefore we 
always have Sk(Ch) — Ch ^ 0. This implies that whenever 
we have a pair (L, C) with L(h) = Sk(Ch) — Ch, we can be 
sure that L is not the zero operator, and we need not worry 
about this requirement any further. 

The analysis in the present case is similar to that carried 
out by Apagodu and Zeilberger [10] , who used it for deriving 
a bound on the order r of L, and similar to our analysis [6] of 
the differential case. The main idea is to follow step by step 
the execution of Zeilberger's algorithm when applied to h. 
This eventually leads to a linear system of equations with 
coefficients in K(n) which must have a solution whenever it 
is underdetermined. The condition of having more variables 
than equations in this linear system is the source of the 
estimate for choices (r, d) that lead to a solution. 

3.1 Zeilberger's Algorithm 

Recall the main steps of Zeilberger's algorithm: for some 
choice of r, it makes an ansatz L — £o+£iS n + - ■ -+£ r Sn with 
undetermined coefficients £q, . . . ,£ r , and then calls Gosper's 
algorithm on L(h). Gosper's algorithm [7] proceeds by writ- 
ing 

S k (L(h)) = Sk(P) Q 
L{h) P S k (R) 

for some polynomials P,Q,R such that gcd(Q, S k (R)) = 1 
for all i G N. It turns out that the undetermined coefficients 
£o,.-.,i r appear linearly in P and not at all in Q or R. 
Next, the algorithm searches for a polynomial solution Y of 
the Gosper equation 

P = Q S k {Y) - RY 

by making an ansatz Y — yo + y\k + yik 1 + ■ ■ ■ + y s k s for 
some suitably chosen degree s, substituting the ansatz into 
the equation, and comparing powers of k on both sides. This 
leads to a linear system in the variables £q, . . . , £ r , yo, . . . , y s 
with coefficients in K(n). Any solution of this system gives 
rise to a telescoper L with the corresponding certificate C = 
RY/ P. If no solution exists, the procedure is repeated with 
a greater value of r. 

For a hypergeometric term h and an operator L — £o + 
£iS n + ■ ■ ■ + £ r S„, we have 



L(h) = Y. £iX 



=i m 

r M 



Aia m jD'b ni 



1 n, 



M 
m — l 



l V r 



p m =iU m m {B, 



m 77i J 



£4**5* (p) II 



X x n y k I 



r(A m )r(B„ 



r([/ m + ru m )Y(V m + rv m ) ' 



where 



Pi, 



■A-m™ Bjn™ {Urn H" VUfi 



(V m + w„ 



\ (r—i)v n 



We can write 



5 fc (L(/i)) _ S fc (P) Q 



L(fc) 



where 



P = £)/ i x\S£(p) P, mi 

i — m— 1 

A/ 

<3 = y TT A^"(v m + rv m — »!)' 



m — 1 
M 



Ji= ]^[ ((7m + ru m - u m Y m B h m 



Depending on the actual values of the coefficients appearing 
in h, this decomposition may or may not satisfy the require- 
ment gcd(<3, Sjfe(-R)) = 1 for all ieK. But even if it does 
not, it only means that we may overlook some solutions, but 
every solution we find still gives rise to a correct telescoper 
and certificate. Since we are interested only in bounding 
the size of the telescopers of h, it is sufficient to study under 
which circumstances the Gosper equation 

P = Q S k (Y) ~ RY 

with the above choice of P, Q, R has a solution. 

3.2 Counting Variables and Equations 

Apagodu and Zeilberger [10] proceed from here by analyz- 
ing the linear system over K(n) resulting from the Gosper 
equation for a suitable choice of the degree of Y. They de- 
rive a bound on r but give no information on the degree d. 
General bounds for the degrees of solutions of linear systems 
with polynomial coefficients could be applied, but they turn 
out to overshoot quite much. In particular, it seems difficult 
to capture the phenomenon that increasing r may allow for 
decreasing d using such general bounds. 

We proceed differently. Instead of a coefficient compari- 
son with respect to powers of k leading to a linear system 
over K(n), we consider a coefficient comparison with respect 
to powers of n and k leading to a linear system over K. This 
requires us to make a choice not only for the degree of Y in k 
but also for the degree of Y in n as well as for the degrees 
of the ii (i = 0, . . . , r) in n. For expressing the number of 
variables and equations in this system, it is helpful to adopt 
the following definition. 

Definition 1. For a proper hypergeometric term h as in (1), 
let 

5 = degp, 

M M 

■d = maxj ^2 ( a "i + bm), ^2 ( Um + w ™)}> 

771—1 m=l 

M 

A = 22 ( Um + w ™)> 



= ~^2 ( am + b 



m (*m 



' »m), 



1V1 M 

maxj ^2 + v' m ), ^22 ( U 'm + O}- 



Note that these parameters are integers which only depend 
on h but not on r or d. Except for /1, they are all nonnega- 
tive. Note also that we have X + fj, > and ■d = A + /i + > 

Lemma 2. Let di := deg n ii (i = 0, . . . , r). Then 
deg P < 8 + Xr + max(cii + in). 

7 — 

Furthermore, deg fe P < 5 + dr. 

Proof. It suffices to observe that 

deg Pi, m < ia m + ib m + (r - i)u m + (r - i)v m 

for all m = 1, . . . , M and all i = 0, . . . , r. For the degree 
with respect to k, observe also that deg fe ii = for alii. ■ 

We have some freedom in choosing the di. The choice 
influences the number of variables in the ansatz 



EE 

i=0 j=0 



ti.jn 1 5L 



as well as the number of equations. We prefer to have many 
variables and few equations. For a fixed target degree d, the 
maximum possible number of variables is (d + l)(r + 1) by 
choosing do = rfi = • • • = d r = d. But this choice also leads 
to many equations. A better balance between number of 
variables and number of equations is obtained by lowering 
some of the di with indices close to zero (if \i is negative) 
or with indices close to r (if fi is positive). Specifically, we 
choose 



di :— d - 



[y + z — r) 

{v-i)+\n\ 



if n > 

if |t < 0. 



See [6, Ex. 11, Ex. 15.5 and the remarks after Thm. 14] 
for a detailed motivation of the corresponding choice in the 
differential case. The support of the ansatz for L looks as in 
the following diagram, where every term n 3 is represented 
by a bullet at position (i,j): 



3=0 



With this choice for the degrees di, the number of variables 
in the ansatz for L is 

r 

Y^{di + l) = (d+ l)(r + 1) - \\n\v{v + 1), 

i=0 

provided that d > \fi\u. The number of resulting equations 
is as follows. 



Lemma 3. If the di are chosen as above, then P contains 
at most 

\(8 + fir + l)(8 + 2d + fir- 2\fi\u + 2) 
terms n l k 3 . 

Proof. If /i > 0, we have 

di + ifj, = d — {u + i — r) fx + ifi < d — + r^i 

for all i = 0, . . . , r. Likewise, when /i < 0, we have 

di + i/i = d — (v — + + i/i < d — 

for all i = 0, . . . , r. Together with Lemma 2, it follows that 

degP < 8 + (A + ^ + )r + d - = 8 + fir + d - v\fi\ 

regardless of the sign of /i. We also have deg fc P < 8 + fir 
from Lemma 2. For the number of terms n % k 3 in P we have 

dog fc P 

(1 + degP-i) = ±(deg fe P+l)(2degP + 2-deg fc P). 

i=0 

Plugging the estimates for deg P and deg fe P into the right 
hand side gives the expression claimed in the Lemma. ■ 

The support of P has a trapezoidal shape which is de- 
termined by the total degree and the degree with respect 
to k: 



dcgP 



The next step is to choose the degrees for Y in n and k. 
This is done in such a way that Q Sk(Y) — RY only contains 
terms which are already expected to occur in P, so that no 
additional equations will appear. 

Lemma 4. Let the di be chosen as before and suppose that 
Y 6 K[n, k] is such that deg Y < degP — v and deg fe Y < 
deg fc P — v. Then P — (Q Sk{Y) — RY) contains at most 

±(<5 + dr + 1) (8 + 2d + dr - 2\fi\v + 2) 
terms n'k 3 . 

Proof. As for Lemma 3, using also maxjdeg Q, deg R} = 
max{deg fe Q, deg k R} — v. m 

Lemma 4 suggests the ansatz 

si s 2 —i 
i=0 i=0 

with si = deg fe P — v and si = deg P — f, which provides us 
with 

| (8 + dr + 1 - v) (8 + 2d + fir - 2\n\v + 2 - v) 

variables. We are now ready to formulate the main result 
of this section. Note that the inequality for d is a consid- 
erably simpler formula than the corresponding result in the 
differential case (Thm. 14 in [6]). 



Theorem 5. Let h be a proper hypergeometric term which 
cannot be written h = qho for some q £ K(n, k) and a hy- 
pergeometric term ho with Sk{ho)/ha — 1. Let 5,X,fJ,,v be 
as in Definition 1, let r >v and 



d> 



($v - l)r + \v(28 + \n\ + 3 - (1 + - 1 



r — v + 1 

Then there exists a telescoper L for h of order r and degree d. 

Proof. A sufficient condition for the existence of a telescoper 
of order r and degree d is that for some particular ansatz, 
the equation 

P = QS k (Y)-RY 

has a nontrivial solution. A sufficient condition for the ex- 
istence of a solution is that the linear system resulting from 
coefficient comparison has more variables than equations. 
For all d in question, we have d > -&v > \\i\v. Therefore, 
with the ansatz described above, we have 

(d + l)(r + l)- H/xKiz+l) 
iP, 

| (8 + -dr + 1 - v) (8 + 2d + dr - 2\n\v + 2 - v) 

variables yij in Y, and 

| (8 + fir + 1) (8 + 2d + -dr - 2\fi\v + 2) 

equations. Solving the inequality 

(d+l)(r + l)- ±|^0 + 1) 

+ | (8 + fir + 1 - v) (8 + 2d + fir - 2\fJ,\p + 2 - v) 

> \{8 + fir + l){8 + 2d + fir- 2\fi\v + 2) 

under the assumption r > v for d gives the claimed degree 
estimate. ■ 

3.3 Examples and Consequences 

Example 6. 1. For 

2 2 £(2n + 3fc) 
h = in +k +1) T(2n-k) 

we have 5 = 2, fi = 2, /j, = 0, v = 4. Theorem 5 
predicts a telescoper of order r and degree d whenever 
r > 4 and 



d> 



7r + 5 



The left figure below shows the curve defined by the 
right hand side (black) together with the region of all 
points (r, d) for which we found telescopers of h with 
order r and degree d by direct calculation (gray). In 
this example, the estimate overshoots by very little only. 

2. The corresponding picture for 

T{2n + k)V{n-k + 2) 
~ T(2n-k)T{n + 2k) 

is shown below on the right. Here, 5 = 0, fi = 3, /i = 0, 
v = 3 and Theorem 5 predicts a telescoper of order r 
and degree d whenever r > 3 and 

, 8r - 1 
d> -. 

r — 2 

In this example, the estimate is less tight. 



d d 




5 10 15 20 5 10 15 20 



The points (r, d) in the portion of the gray region which is 
below the black curve represent telescopers where the corre- 
sponding linear system resulting from the ansatz considered 
in our proof is overdetermined but, for some strange rea- 
son, nevertheless nontrivially solvable. The small portions 
of white space which lie above the curves are not in contra- 
diction with our theorem because they do not contain any 
points with integer coordinates. (The theorem says that ev- 
ery point (r, d) € 71? above the curve belongs to the gray 
region.) 

Theorem 5 supplements the bound given in [10] on the or- 
der of telescopers for a hypergeometric term by an estimate 
for the degree that these operators may have. In addition, it 
provides lower degree bounds for higher orders and admits 
a bound on the least possible degree for a telescoper. 

Corollary 7. With the notation of Theorem 5, h admits a 
telescoper of order r = v and degree 

as well as a telescoper of order 

r= [^(1 + 2(5 + 2(1/- 
and degree d — -&v . 

Proof. Immediate by checking that the two choices for r and 
d satisfy the conditions stated in Theorem 5. ■ 

An accurate prediction for the degrees of the telescopers 
can also be used for improving the efficiency of creative tele- 
scoping algorithms. Although most implementations today 
compute the telescoper with minimum order, it may be less 
costly to compute a telescoper of slightly higher order. If 
we know in advance the degrees d of the telescopers for ev- 
ery order r, we can select before the computation the order r 
which minimizes the computational cost. Of course, the cost 
depends on the algorithm which is used. It is not necessary 
(and not advisable) to follow the steps in the derivation of 
Theorem 5 and do a coefficient comparison over K. Instead, 
one should follow the common practice [8] of comparing co- 
efficients only with respect to powers of k and solve a linear 
system over K(n). For nonminimal choices of r, this sys- 
tem will have a nullspace of dimension greater than one, of 
which we need not compute a complete basis, but only a 
single vector with components of low degree. There are al- 
gorithms known for computing such a vector using 0(m 3 t) 
field operations when the system has at most m variables 
and equations and the solution has degree at most t [3, 13, 
5]. In the situation at hand, we have m — (r+l) + (<5+i?r+l) 
variables and a solution of degree t = S+'&r+d— {\ij,\+1)v+1. 



Therefore, in order to compute a telescoper and its certifi- 
cate most efficiently, we should minimize the cost function 

C(r,d) ~ (($ + l)r + S + 2) 3 (S + i}r + d-(\Li\ + l)v + l). 

According to the following theorem, for asymptotically large 
input it is significantly better to choose r slightly larger than 
the minimal possible value. 

Theorem 8. Let h and A, fi, v be as in Theorem 5, r > 
max{i?,j/}, and suppose that k £ 1R is a constant such that 
degree t solutions of a linear system with m variables and 
at most m equations over K(n) can be computed with nm 3 t 
operations in K. Then: 

1. A telescoper of order r = r along with a corresponding 
certificate can be computed using 

«r 9 + i(7-|At|)«;T 8 +0(r 7 ) 

operations in K. 

2. If a > 1 is some constant and r is chosen such that 
r — or + O(l), then a telescoper of order r and a 
corresponding certificate can be computed using 



operations in K. 

In particular, a telescoper for h and a corresponding certifi- 
cate can be computed in polynomial time. 

Proof. According to Theorem 5, for every r > r there exists 
a telescoper of order r and degree d for any 

d>f(r) := fr'-fr + Ofr'). 

r — t — 1 

By assumption, such a telescoper can be computed using no 
more than 

C(r, d) = k((t + l)r + S + 2) 3 (S + rr + d- (\fj,\ + l)r + l) 

operations in K. The claim now follows from the asymptotic 
expansions of C(r, /(r) + l) and C(ar, /(ar) + l) for r — > oo, 
respectively. ■ 

The leading coefficient in part 2 is minimized for a — 5/4. 
This suggests that when § and v are large and approximately 
equal, the order of the cheapest telescoper is about 20% 
larger than the minimal expected order. 

4. THE RATIONAL CASE 

We now turn to the case where h can be written ash = qho 
for some hypergeometric term ho with Sk(ho)/ho = 1. By 
the following transformation, we may assume without loss 
of generality ho = 1. 

Lemma 9. Let h be a hypergeometric term and suppose 
that h = qho for some q £ K(n, fc) and a hypergeometric 
term ho with Sk(ho)/ho = 1. Let a,b £ K[n, k] be such that 
S n (ho)/ho = a/6. Let L be a telescoper for q of order r and 
degree d. Then there exists a telescoper for h of order r and 
degree at most d + r max{deg n a, deg n b}. 



Proof. Write L = £ +£iS n + ■ ■ ■ +l r S r n and let C G K(n, k) 
be a certificate for L and q, so L(g) = Sk(Cq) — Cq. For 
i = 0, . . . , r, let 




and L := i + hs n + ■■■+ l r S r n . Then 

L{qh ) = L{q)h = {S k (Cq) - Cq)h = S k (Cqh ) - Cqh . 
Because of 

(0\ = S k (S n (h )) = Sn(g fc Qo)) = ■S'n(feo) = O 

the operator L is free of k. Thus, after clearing denomina- 
tors, L is a telescoper for h with coefficients of degree at 
most d + r max{deg n a, deg n b}. m 

From now on, we assume that h is at the same time 
a proper hypergeometric term and a rational function, or 
equivalently, that h is a rational function whose denomina- 
tor factors into integer-linear factors. Le [9] gives a precise 
description of the structure of telescopers in this case, and 
he proposes an algorithm different from Zeilberger's for com- 
puting them. Our degree estimate is derived following the 
steps of his algorithm, so we start by briefly summarizing 
the main steps of Le's approach. 

4.1 Le's Algorithm 

Given a rational proper hypergeometric term h, Le's al- 
gorithm computes a telescoper L for h as follows. 

1. Compute g G K(n, k) and polynomials p,q G K[n, k] 
with gcd(c/, Sl(q)) = 1 for alH € % \ {0} such that 

P 

h = Sk(g) - g + -■ 

Then an operator L is a telescoper for h if and only 
if L is a telescoper for 2. Abramov [1] and Paule [11] 
explain how to compute such a decomposition. 

2. Compute a polynomial u G K[n], operators Vi, . . . , V s 
in K[n][S„], and rational functions fi,...,f a of the 
form fi = (can + a'jk + a")~ e ' (i = 1, . . . , s) such that 

* i=0 

Such data always exists according to Lemma 5 in [9] 
in combination with the assumption gcd(g, S k (q)) = 1 
(i G Z5 \ {0}). It can be further assumed that the fi 
are chosen such that a\ > 0, e t > 0, gcd(ai,cii) = 1 for 
all i, and 

for all i ^ j with ei = ej. 

3. For i = 1, . . . , s, compute an operator Li G K(n)[S n ] 

such that Sn z — 1 is a right divisor of Li(-Vi). It follows 
from Le's Lemma 4 that the operators Li with this 
property are precisely the telescopers of the rational 
functions Vi(fi). 

4. Compute a common left multiple L G K[n][Sn] of the 
operators Li, . . . , L s . Then L is a telescoper for h. 



The main part of the computational work happens in the 
last two steps. It therefore appears sensible to assume in 
the following degree analysis that we already know the data 
u,V±, . . . ,V 3 , fi, . . . , fs computed in step 2, and to express 
the degree bounds in terms of their degrees and coefficients 
rather than in terms of the degrees of numerator and de- 
nominator of h, say. 

4.2 Counting Variables and Equations 

Also in the present case, the degree estimate is obtained 
by balancing the number of variables and equations of a 
certain linear system over K. The linear system we consider 
originates from a particular way of executing steps 3 and 4 
of the algorithm outlined above. 

Theorem 10. Let u G K[n] and let Vi,...,V 3 G K[n][S„] 
be operators of degree Si (i = l...,s). Let ft — {ain + 
a[k + a")~ ei for some a" G K, ai,a'i G % with a[ > and 
gcd(ai,a'i) = 1, ei > 0, suppose 

Vet; a' 3 J V< arJ 

for all i 7^ j with e, = ej. Let h = — Y2t=i Vi(fi)- Then for 
every r > Xli=i a i an d ev ^ry 

3 

-r—1 + y^^a'jSi 

d > ^ h deg n u 

r + 1 - a'i 

there exists a telescoper L for h of order r and degree d. 

Proof. According to Le's algorithm, it suffices to find some 
L G KfnJfS'n] and operators Ri G TK.(n)[S„] with the prop- 
erty that L(iVi) = RiiSn* - 1) for all i. 

Denote by pi the order of Vi. Writing d := d — deg n u, we 
make an ansatz L — Lu with 

r d 
i=0 j = 

so that L has degree d and L^Vi = LVi (i = l,...,s). 
It thus remains to construct operators Ri G K[n][S n ] with 

LVi — Ri(S„ i — 1). Since LVi has order r + pi and degree 
d + Si, we consider ansatzes for the Ri of order r + pi — a'i 

and degree d + 8i, respectively, because Sn z — 1 has order a'i 
and degree 0. Then we have altogether 

s 

(r + l)(d+ 1) + J2( r + P>-^ + + & + 1) 
i=i 

variables in L and the and comparing coefficients with 
respect to n and S n in all the required identities LVi = 
Ri(S„ l — 1) leads to a linear system with 

s 

Y,( r + P> + + & + 1) 

i=i 

equations. This system will have a nontrivial solution when- 
ever the number of variables exceeds the number of equa- 



tions. Under the assumption r > a i> the inequality 

(r + 1) + Y^{r + pi -a'i + l)(d + & + 1) 

1=1 

>^(r+pi + l)(d + 5» + l) 

i=l 

is equivalent to 



d> 



-r - 1 + gift 



This completes the proof. ■ 

4.3 Examples and Consequences 

Example 11. 1. The rational function 
(2n - 3fc)(3n - 2fc) 2 



ft = 



(n + k + 2)(n + 2fc + l)(2n + fc + l)(3n + k + 1) 

can fee written in the form h = — Y^t=i Vi(fi) where 
u = (n-l)n(n + 3)(2n-l)(3n+l)(5n+l), t/ie /» are 
smc/i f/iat = a' 2 = a-j = 1, a' 4 = 2, and i/ie Vi are 
such that Si = ■ ■ • = 84 = 6. Therefore, Theorem 10 
predicts a telescoper of order r and degree d whenever 
r > 5 and 

, 29 - r „ 

d > +6. 

r — 4 

TTiis curve together with the region of all points (r, d) 
for which a telescoper of order r and degree d exists is 
shown in the left figure below. 

2. The corresponding picture for the rational function 



h 



(n-fc + l) 2 (2n-3fc + 5) 



(n + k + 3)(n + k + 5)(n + 2k + l)(2n + k + l) 2 



is shown in the figure below on the right. This input 
can be written in the form 



1 4 

h = S k (g)-g+-J2v i (f i , 



with _ 10(n+3) 2 (, 1+ 4)(2n+2fc+7) _ , -,-,2 / _ 

4) 2 (n - 2) 2 (n + 5)(n + 9), the fi such that a[ = a' 2 = 
a'j = 1, a' 4 = 2, and the Vi such that Si = 8, S 2 — 5$ = 
<$4 = 7. According to Theorem 10, we therefore expect 
a telescoper for h of order r and degree d whenever 
r > 5 and 

, 35 - r a 
d > +8. 

r — 4 

In this example, the estimate is not as tight as in the 
previous one. 



Again, it is an easy matter to specialize the general degree 
bound to a degree estimate for a low order telescoper, or to 
an order estimate for a low degree telescoper. 

Corollary 12. With the notation of Theorem 10, h admits 
a telescoper of order r — YjI=i a 'i ana " degree d — deg n u + 
Si=i(& — as we ^ as a telescoper of order r — Y2t=i a ^ 
and degree d = deg n it. 

Proof. Clear by checking that the proposed choices for r 
and d are consistent with the bounds in Theorem 10. ■ 

Also like in the non-rational case, the bounds for the de- 
grees of the telescopers can be used for deriving bounds on 
the computational cost for computing them. In the present 
situation, let us assume for simplicity that the cost of steps 
1 and 2 of Le's algorithm is negligible, or equivalently, that 
the input h is of the form i Yli=o ^(/O with Vi € TK[n][S n ]- 
We shall analyze the algorithm which carries out steps 3 
and 4 of Section 4.1 in one stroke by making an ansatz over 
K(n) for an operator L — £o + £iS n + ■ ■ ■ + irS^, comput- 
ing the right reminders of L^Vi with respect to S„' — 1 and 
equating their coefficients to zero. We assume, as before, 
that the resulting linear system is solved using an algorithm 
whose runtime is linear in the output degree and cubic in the 
matrix size. Then the algorithm requires 0(r 3 d) operations 
in K. 

Theorem 13. Let u € K[n], Vi,...,V s £ K[n][S n ], and 
fi, . . ■ , fs £ K(n, k) be as in Theorem 10 and consider h — 
u Si=i Vi(fi)- Suppose that k £ ]R is a constant such that 
degree t solutions of a linear system with m variables and 
at most m equations over K(n) can be computed with nm 3 t 
operations in K. Assume 5i = ■ ■ ■ = 8 S —: S > and 
ai = 02 = • ■ ■ = a s =: a are fixed. Then: 

1. A telescoper of order r = as can be computed using 

a' 4 (5-l) K s 4 + 0(s 3 ) 
operations in K. 

2. If a > 1 is some constant and r is chosen such that 
r = aa's + 0(1) then a telescoper of order r can be 
computed using 

3 

-a 3 (5 — 1 + (a — 1) deg n ii)ks 3 + Ofs 2 ) 

a — 1 

operations in K. 

In particular, a telescoper for h can be computed in polyno- 
mial time. 





5 10 15 20 05 10 15 20 



Proof. According to the first estimate stated in Theorem 10, 
for every r > a's there exists a telescoper of order r and 
degree d for any 



d > f(r) := 



z'<5- 



1 



r + 1 



deg n u. 



By assumption, such a telescoper can be computed using no 
more than C(r, d) := nr 3 d operations in K. The claim now 
follows from the asymptotic expansions of C(a's, f(a's) + 1) 
and C(aa's, f(aa's) + 1) for s — t 00, respectively. ■ 



When deg n u — 0, the leading coefficient in part 2 is min- 
imized for a = 3/2. This suggests that when s is large and 
all the Si, and a\ are approximately equal, the order of the 
cheapest operator exceeds the minimal expected order by 
around 50%. 

It must not be concluded from a literal comparison of the 
exponents in Theorems 5 and 10 that Le's algorithm is faster 
than Zeilberger's, because r in Theorem 5 and s in Theo- 
rem 10 measure the size of the input differently. Neverthe- 
less, it is plausible to expect that Le's algorithm is faster, 
because it finds the telescopers without also computing a 
(potentially big) corresponding certificate. Our main point 
here is not a comparison of the two approaches, but rather 
the observation that both of them admit a degree analysis 
which fits to the general paradigm that increasing the order 
can cause a degree drop which is significant enough to leave 
a trace in the computational complexity. 

It can also be argued that the situations considered in 
Theorems 8 and 13 are chosen somewhat arbitrarily (i? and 
v growing while pi remains fixed; resp. s growing while all 
the Si and a' t remain fixed). Indeed, it would be wrong to 
take these theorems as an advice which telescopers are most 
easily computed for a particular input at hand. Instead, in 
order to speed up an actual implementation, one should let 
the program calculate the optimal choice for r from the de- 
gree estimates given Theorems 5 and 10 with the particular 
parameters of the input. 

Unfortunately, we are not able to illustrate the speedup 
obtained in this way by an actual runtime comparison for a 
concrete example, because for examples which can be han- 
dled on currently available hardware, the computational cost 
turns out to be minimized for the least order operator. But 
already for examples which are only slightly beyond the ca- 
pacity of current machines, the degree predictions in The- 
orems 5 and 10 indicate that computing the telescoper of 
order one more than minimal will start to give an advan- 
tage. We therefore expect that the results presented in this 
paper will contribute to the improvement of creative tele- 
scoping implementations in the very near future. 
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